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Abstract. We consider the KPZ equation in one space dimension with narrow 
wedge initial condition, h(x,t = 0) = —\x\/5, i5 < 1, Based on previous results 
for the weakly asymmetric simple exclusion process with step initial conditions, we 
obtain a determinantal formula for the one-point distribution of the solution h(x,t) 
valid for any x and t > 0. The corresponding distribution function converges in 
the long time limit, t — > oo, to the Tracy- Widom distribution. The first order 
correction is a shift of order t~ 1//3 . We provide numerical computations based on the 
exact formula. 



1 Introduction 



In a seminal paper [l|| Kardar, Parisi, and Zhang (KPZ) proposed a continuum 
equation for surface growth, which in the particular case of one space-dimension 
reads 

Here h(x, t) is the height profile at time t, t > 0, and x£l. A is the strength of the 



nonlinearity. (1.1 ) is invariant under A ~^ — A and h —h. For later convenience we 
choose A > 0. v > is the coefficient of diffusive relaxation, rj normalized Gaussian 
white noise with covariance (rj(x,t)r}(x',t')) = 8{x — x')5(t — t'), and \[D the noise 



strength. (1.1) is the KPZ equation in one space dimension. 

Typical realizations of the white noise rj(x, t) are far from regular functions. 
The solution h(x,t) partially inherits this roughness of the noise and therefore the 
square (dh/dx) 2 is ill-defined, in general. Nevertheless meaningful solutions can be 
constructed by suitable approximation schemes. They will be explained in detail 
below. The most direct one can be easily stated. One smoothens r\ to r\ K as 

T} K (x, t) = J <\x'kip(k(x - x'))r](x', t) = cp K * r](x, t) (1.2) 

with ip some smooth, localized, and normalized smearing function. Then r] K — > rj as 



k — y oo. (1.1) with noise r\ K has well-defined solutions, denoted by h K (x,t). They 
move with a uniform background velocity v K along the /i-direction. v K — > oo as 
k oo, but h K (x, t) — v K t has a limit. Since v K sets merely the choice of a reference 
frame, the claim is that under this limit procedure the fluctuation properties remain 
intact. 

While meaningful solutions are thus ensured, very little is known about their 
properties. To make dh(x, t)/dx stationary, one has to start the KPZ equation with 
two-sided Brownian motion. With this input, one argues that the height fluctuations 
will grow as t 1//3 , while the transverse correlation length increases as 

£2/3 very 

recently it has been proved that the variance of the stationary two-point function 
increases as t 4//3 by providing suitable upper and lower bounds 2 . 



In our investigation we choose as initial data the narrow wedge 

h(x, 0) = -\x\/6 (1.3) 

with 5 < 1. This may look artificial. However, for short times the nonlinearity 
dominates and, ignoring the other terms in the equation, h spreads very rapidly 
into the parabolic profile 

h(xt) = h 2/2Xt f° r l x l — 2 At/5 , 
V ' ; \-\x\/S for |x| > 2M/6, ' 

which should be viewed as the top part of a growing droplet. Physically one thereby 
covers the case of macroscopically curved height profiles. 
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We will provide a determinantal formula for the one-point distribution of h(x, t) 
at prescribed x and valid for all t > 0, which is exact in the properly normalized 



limit 5—7-0. The formula will be given in Equation (4.17). Numerical plots based 
on this formula are provided in Section 6. 

The exact solution is constructed from the corner growth model, which we intro- 
duce first. The corner growth is the stochastic evolution of a height function h(j, t), 
j G Z, t > 0, which takes integer values constrained by \h(j + l,t) — h(j,t)\ = 1. 
The initial condition is the wedge h(j,0) = — In a short time interval dt, inde- 
pendently for each height, at a local minimum the height function increases by 2 
with probability pdt and stays put with probability 1 — pdt. Correspondingly, at a 
local maximum it decreases by 2 with probability qdt and stays put with probability 
1 — qdt. These rules respect the constraint. We set q + p = 1 and < p < q < 1, so 
the height will decrease on average. It has been proved in |4j|5), see also [6] for the 
case p — 0, that 

h(0, t) = -\{q-p)t + 2- 1 / 3 ((g - p)tY^ TW (1.5) 

for large t. 1/3 is the scaling exponent predicted already on the basis of the KPZ 
equation ( |1.1 ). £ TW is a random amplitude which is Tracy- Widom distributed 



i.e. 



P(£tw <s) = Ftw(s) = det(l - P s K Ai P s ) (1.6) 
with K Ai the Airy kernel, 

POO 

K Ai (x,y) = dAAi(x + A)Ai(|/ + A), (1.7) 
Jo 

Ai the Airy function, and with P s the projection onto the interval [s, oo). P S K A \P S 
is a trace class operator on L 2 (IR). Hence the Fredholm determinant (1.6) is well- 
defined on L 2 (R). 

As will be discussed in detail, the exact KPZ solution is obtained through the 
corner growth model in the limit of small asymmetry, more precisely one requires 
q—p 1, space of order (q — p)~ 2 , and time of order (g — p)~ 4 . For large t the exact 



solution attains the Tracy- Widom statistics as in (1.5). The solution also predicts 
that the first order correction to Tracy- Widom is a shift of order t~ 1//3 . Such quali- 
tative features remain valid over a much wider range of q — p than anticipated by the 
method of construction |3]. For the corner growth model the maximal asymmetry 
is q = 1. We expect that upon decreasing q towards 1/2 finer and finer properties of 
the exact solution are revealed. At present such claims can be tested only through 
Monte Carlo simulations. For example, at q — 1 for t = 10 3 MC steps the average 
height is —512 and its distribution has an effective range of 30 steps each of size 2. 
Even at this discrete level the Tracy- Widom distribution is an accurate approxima- 
tion provided one includes a shift of 2.9 to the right relative to the approximation of 
(1.5). For smaller values of q, approximately at q = 0.78, the prefactor of the t -1 / 3 



shift changes sign, which again can be understood from the KPZ solution. To leading 
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order in t the solution is shifted by 7 t ~ 1 (— 0.577+2 log(g— p)), 7$ = 2~ 1 / 3 (g— p) 4//3 t 1//3 , 
which is indeed negative. To study finer details would require longer Monte Carlo 
runs, which have not been carried out yet. Available is the t -1 / 3 correction of the 
Tracy- Widom distribution for the polynuclear growth model, as based on an exact 
recursion formula [7]. For large t a definite limit is approached. The correction has 
a single node. Thus it agrees qualitatively with the t -1 / 3 correction of our solution, 
which is proportional to F^ w (s), but the precise shape of the correction is different. 
Apparently the first order correction to Tracy- Widom depends on the details of the 
model. 



2 The Cole-Hopf solution 

We start from the mollified KPZ equation 

dt K = ^ X{ dx K)2 + U d^ K + ^ ' (2 - 1} 



where the Gaussian noise r/ K is smeared as in <\1.2 ) and hence has the covariance 



(Vk(x, t)r} K (x', £')> = <p K * <p K (x - x')5(t - if) . (2.2) 
The mollifier ip K satisfies (p K (x) = kip(kx), (p > 0, y? of rapid decay, <p(x) = ^p(-x) 



and j ip(x)dx = 1. The average over space-time white noise is denoted by (•). (2.1) 
becomes linear through the Cole-Hopf transformation 

Z K (x, t) = exp[(X/2v)h K (x, t)} , (2.3) 

in which case Z K satisfies 

-Z K = Vq^Z k + (\VD/2v) VK Z K . (2.4) 



The solution to (2.4) can be written in terms of the Feynman-Kac formula by intro- 
ducing an auxiliary standard Brownian motion b(t) with variance t. Denoting K x 
its expectation when starting at x, one has 

Z K (x,t)=E x (exp[(\>/D/2v) J r) K (b(2vs),t - s)ds]Z(b(2ut), 0)) (2.5) 

for some deterministic initial condition Z(-,0). 

In our context, the initial condition is a narrow wedge of the form (1.4) with a 
global off-set. More precisely we assume h(x,0) = (2u/X)(— — log 25), S <C 1. 
Then the Cole-Hopf transform Z(x, 0) = {25)~ 1 exp[— \x\/S] which is sharply peaked 
at x = and, as an idealization, we will set 

Z(x,0) = 5(x). (2.6) 
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Note that the average 



(Z K (x,t)) = exp [\{\VD/2ufip R * tp K (Q)t]m x (Z(b(2ut), 0)) . (2.7) 

Since ip K * ip K (0) — k, it diverges as exp[c/t] with a suitable constant c > 0. Thus to 
have a limit as k — >■ oo, one better divides the right hand side of (2.5) by (Z K (x,t)), 
which corresponds to the Wick ordering of the exponential, see for example |10|. 
Hence we define 

:Z K {x,t): = Z K (x,t)exp[-l(\VD/2v) 2 <p K *<p K {0)t] . (2.8) 



It is proved in [9] that, for continuous Z(-,0) with an explicit decay condition for 
\x\ — > oo, the limit 

lim :Z K (x,t): = Z(x,t) (2.9) 

exists. Z(x,t) is a well-defined stochastic process. With probability one it has 
sample paths which are continuous in both x, t. In addition Z(x, t) > provided 
Z(x, 0) > and Z(-,0) ^ 0. With this input we define 

h(x,t) = -^-logZ(x,t) (2.10) 
A 

as the solution of the KPZ equation with 5-correlated noise. For later purposes it is 
convenient to introduce, for arbitrary a > 0, 

ah a {x,t) = logE x (:exp [a J dsrj(b(s),t - s)]: 5(6(t))) . (2.11) 

Clearly h a is the solution to the KPZ equation with D = 1 and a = X/2u. By 
construction the first exponential moment is centered as 

(e Q M*,*)) = (27rt)- 1 / 2 e- (:E2/2 ' ) . (2.12) 

From the point of view of the height function, the subtraction corresponds to 

h K (x, t) - \(\D/2u)^ K * if K {0)t . (2.13) 

Hence the uniform velocity, alluded to in the introduction, is given by 

v K = ±(\D/2v)<p K *<p K {0). (2.14) 



(2.5 ) can be viewed also as a directed polymer subject to a space-time white noise 
random potential. With the initial condition (2.6) the polymer b(s) starts at x and 
ends at at time 2vt, which corresponds to a point-to-point directed polymer. The 
potential of the random environment is somewhat singular, hence the subtraction. 
Z(x,t) is the random partition function of the directed polymer at the "inverse 
temperature" \\[Dj2v. From the perspective of disordered systems, of interest is 
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the random free energy (2i>/ A) log Z(x, t), which is nothing else than the KPZ height 
h(x, t). 

The observation leading to (2.7) extends to higher moments. We consider the 
n-th moment of :Z K (t): and for this purpose introduce n independent standard 
Brownian motions bj(t), j — 1, . . . , n, as replicas of b(t). Then 

n „i 

((:Z K (x,t):) n )=E x x . . . x E x (( exp [(Av^D/2z/) V / VkM^us), t - s)ds 

V 3=1 J ° 

n 

-\n(\\fDl2v) 2 ip K * <p K (0)t] > J] Z(b J (2ut), 0) 

3=1 

n ,t 

EJexp [\{\y/D/2u) 2 ^ / 

i,i=lM.i Jo 



E x x ... x 



n 

x JJZ(6 i (2i4),0) 



(2.15) 



which is the path integral for n quantum particles with pair potential —ip K *ip K . Thus 
the Wick ordering of the exponential in (2.9 ) is identical to the normal ordering of the 
pair interaction. (p K (x) — > 6(x) as k -> oo and (2.15) converges to the path integral 
of a system of n quantum particles on the line interacting through an attractive 
(^-potential. The corresponding hamiltonian reads 



n 1 d 2 1 n 

3=1 J hj=l^¥=j 



(2.16) 



where 



a = {2u)-' i,2 \D l/2 . (2.17) 

Since Q~ t2vH ^ has a continuous kernel, it is a meaningful limit to set Z(x, 0) = 5(x) 
in (2.15). Let us denote the n-particle state YYj=i $( x ~ x j) ^y Then 

(Z{x,t) n ) = (0\e- t2 " Hn \x) (2.18) 



with inner product in L 



For later purposes we record that, for the initial condition (2.6), 

1 



(:Z K (x,t/2v):) 



'2nd 



exp[-x 2 /2t] 



(2.19) 



and hence also in the limit k — >■ oo, in agreement with (2.18) for n = 1. 

Very recently, Balazs, Quastel, and Seppalainen [2] studied the KPZ equation 
with initial condition Z(x,0) = exp[(X/2u)b(x)], where b(x) is two-sided Brownian 
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motion. They consider the solution as defined in (2.10) and prove, amongst other 
related results, upper and lower bounds of the form 



c_t 2/3 < Var (h(x, £)) < c + t 2/3 . (2.20) 

We cannot treat their case, because we miss the yet to be accomplished extension 
of the contour integration formula in |4j to initial conditions given by the Bernoulli 
| measure, see M for recent progress. 



3 The WASEP limit 



Bertini and Giacomin [9] (for short BG) discovered a rather different, and in a way 
physically better motivated, limiting procedure for the construction of solutions to 
the KPZ equation. They consider the weakly asymmetric version of the simple 
exclusion process (WASEP), which is equivalent to the corner growth model alluded 
to already in the introduction. In the height picture this corresponds to the following 
stochastic evolution: A height function h : Z — > Z has to satisfy the constraint 

\h(j + 1) - h(j)\ = 1 (3.1) 

for all j G Z. Independently, each height variable waits a unit exponentially dis- 
tributed time. Then h(j) ~» h(j) — 2 with probability q and h(j) ~* h(j) + 2 with 



probability p, p + q — 1. Either transition is suppressed in case (3.1) is violated. 



The height differences are governed by the partially asymmetric simple exclusion 



process (PASEP). Here, because of (3.1), there is at most one particle per site and, 
under the exclusion constraint, particles jump with rate p to the right and rate q to 
the left. We choose q > p. Thus h decreases to — oo. 

Weak asymmetry means that q — p is small. The precise asymptotics can be 
studied in a scaling limit with dimensionless scale parameter e, e <C 1. We choose a 
diffusive space-time scaling, i.e. time as e~ 2 t, t = Oil), space as [e -1 ^], iEK, ['J 
denoting integer part. The choice q — p = 0(e) is the widely studied WASEP, which 
has the deterministic KPZ equation (D = 0) as limit and Gaussian fluctuations 



relative to that profile 11 . In contrast to the conventional choice BG consider 



the crossover scale with an asymmetry equal to y^e, see 12 for a more complete 



discussion. For purpose of comparison with the KPZ equation we set more generally 

q-p = /3^, 0>O. (3.2) 

For the initial height profile BG assume 

h(j,0)= Le- 1/2 0fc7')J, (3.3) 

where is a smooth function with some decay condition at infinity. Here |_"J means 



integer part subject to (3.1 ). We denote by h £ (j, t) the WASEP time-evolved height 
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profile, the superscript e reminding that the rates are e-dependent. BG prove the 
limit 

^Ve(/i e (Le -1 a;J,e- 2 t) + \te^ /2 ) = h(x,t) + ±t , (3.4) 

where h(x, t) is the solution of the KPZ equation, as defined in the previous section, 
with initial condition (p and parameters X — 1, u — ^, D — 1. 

Note that h(j + l, 0)—h(j, 0) = 0(y/e) and the particle density hence deviates by 
order y/e from |. At that density the average time integrated particl e cur rent equals 



— fc~ 3 / 2 /4 to leading order, which is the subtraction on the left of (3.4). (There is 
a factor 2 when switching from integrated currents to heights.) The subtraction on 
the right is more subtle. As will be discussed in more detail below, it originates 
from the normalization condition 

( e M*,0) = f (2vrt)- 1 / 2 exp[-(x - x') 2 /2t]e^dx' . (3.5) 



To implement for the WASEP height the narrow wedge initial condition from 
Section 2, the natural choice would be density \ — S^yfe to the left and density 
| + d^yfe to the right of the origin. Unfortunately, for such initial conditions there 
are no manageable formulas available and we have to take resort to initial conditions 



which no longer are covered by (3.4) and the results in |9j. 

As standing assumption, the WASEP has asymmetry /3y/e and an initial config- 
uration such that all sites from 1 to oo are filled and all other sites are empty. The 
corresponding initial height profile is h(j,0) = After some short initial time 

span, the typical height profile becomes 

h°(ie-H)<*t -Wl^-\$te-^ forb-|</3te-^, 



compare with (1.4). If one restricts j to the interval [— e 1 ,e 1 ], then there one has 



a parabolic profile, which scales exactly as assumed in (3.3) and thus supports that 



(3.4) remains valid for the initial sharp wedge. In 12 we obtained the one-point 
height distribution for the WASEP with 0-1 initial condition in the limit e — > 0. 
Thus the plan is to use this result in the construction of the one-point distribution 
of the KPZ equation with narrow wedge initial condition. 

Let us denote by T)j(t) = 0, 1 the occupation variables of the WASEP at site j 
and time t. Since there is a random index j such that T)j(t) = for j < j , one can 
define the WASEP height by 

3 

h%j,t) = -2 ^ vS)+j. (3.7) 

e=-oo 

Let JT £ (j,t) be the particle current across the bond + 1) integrated up to time 
t. Then 

h e (j,t)=2J e (j,t)-\j\. (3.8) 
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In 12 we studied the distribution of J e at j = [y(3te~ 3 ^ 2 + e~ l x\ and time s~ 2 t. To 
identify the height statistics of the KPZ equation it suffices to take y = 0, in which 

case 



limP(^/i £ (L£ _1 xJ,£- 2 + \(5He- 1 + (x 2 /2t) - \og(2^) < 7i s) 

= F t (s) , (3.9) 

with 

7f = 2" 1 / 3 (/3 4 t) 1 / 3 . (3.10) 
F t (s) is a t-dependent family of distribution functions. To define them, let us 



introduce the Airy kernel K Ai (x,y) of (1.7), the kernel 

Bt(x,y) = K Ai (x,y) 

POO 

+ / dA(e 7 ' A - I)" 1 (Ai(x + X)Ai(y + A) - Ai(x - X)Ai(y - A)) , (3.11) 
Jo 

and the unnormalized projection P A i with kernel 

P Ai (x,y)=Ai(x)Ai(y). (3.12) 

Then 

/oo 
dwexp[-e 7t(s - M) ] 
-oo 

x ( det(l - P u (B t - P Ai )P u ) - det(l - P u B t P u )) . (3.13) 
Here P u projects onto [u, oo). P u B t P u and P u P Ai P u are trace class operators in 



L 2 (IR). Hence the Fredholm determinants in (3.13) are well-defined. 

Note that the probability density 4-F t (s) is the convolution of the Gumbel prob- 
ability density 7te 7tS exp[— e 7 * s ] and the difference of determinants, 

g t (u) = det(l - P u (B t - P Ai )P u ) - det(l - P u B t P u ) . (3.14) 

For t — t- oo the Gumbel density converges to 5(s) and gt(u) to F^ w (u), see (l~2] . 
Hence 

lim F t (s) = Ftw(s) (3.15) 

i— >-oo 



with -Ftw the Tracy-Widom distribution function (1.6) 



4 Centering and exact solution 

Through the limit to 5-correlated noise we adjusted the first moment of the partition 
function as 

{Z(x,t/2u)) = -^=exp[-^] 1 (4.1) 



L 2t- 
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which sets an important constraint for the WASEP limit. To exploit it, one notes 
that a particular exponential moment of the WASEP also satisfies a closed linear 
equation. This property was first proved by Gartner 13 , 14 , see also [9] and the 
related contribution 15 . For the PASEP with fixed p, q and arbitrary initial con- 
ditions, let us denote by E t the expectation with respect to the height statistics at 
time t and let us set 

/(j,t)=E t (e^). (4.2) 

"d is adjusted such that 



e~ 2 » = V - . (4.3) 
Q 



Then / is the solution of 



Jt M 4) = =s* G'G + 1 ' t) + M " M) " f(j - 4)) + < Jhi " mM • 

/(j,0)=e-*W. (4.4) 



(4.4) can be deduced by working out the time derivative as 

|/(j,t)=E t (Le^)) (4.5) 

with L the generator of the Markov jump process for the heights. To satisfy the 
constraint (4.1) one has to consider the height shifted by log / as 

$h e ( [e~ x x\ , e~H) - log /( [e~ 1 x\ , £~ 2 t) - (x 2 /2t) - log V2^i . (4.6) 

/ can be written as Fourier integral, 

f^xle-H) 

i r 

= — J dfcexp [(cos k — l)(cosh$) -1 £ -2 £ + ike~ 1 x 

+ \ogf{k,0)+e- 2 t((cosh$y 1 - 1)] (4.7) 

with the initial condition 

f(k,0) = (1 — e _2,? )(l — 2e~® cos k + e _2l? ) _1 . (4.8) 
log/(|_£ _1 2;J ,s~ 2 t) has to be computed including terms of order 1. One has 

= pe 1 ' 2 + |/3 3 e 3 / 2 + 0{e h ' 2 ) , (4.9) 

(coshtf)- 1 = 1 - \fe - \^e 2 + 0{e 3 ) . (4.10) 

The initial condition (4.8) is approximated by 2//3y/e. Then the Fourier integral is 
the transition probability of a simple random walk with rate (cosh0) _1 starting at 0, 
which at the given scale can be approximated by a normalized Gaussian. Therefore 

-log/(le-y|,e- 2 t) 



pHe- 1 + (x 2 /2t) + \pH + log(V2^/e) + log(/3v^/2) + 0(y/i) . (4.11) 
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Inserting in (4.6) and taking the deterministic approximation to h £ for the second 
order in the expansion of one arrives at 

dh £ ( [rtj , e'H) - log /( [e~ 1 x\ , e~ 2 t) - (x 2 /2t) - log V2^i 

= f3^eh £ ( [e-^J , £- 2 t) + l^te" 1 - ^/3 4 t - log(2v^//3) + 0( ^i) • (4.12) 



Note that the e-dependent part of the centering is precisely as in (3.9). 

For the WASEP at the reference point y = the average particle density p equals 
1/2 and the noise strength equals p(l — p) = 1/4. Including the factor 2 for the 
conversion between density and height differences, the corresponding parameters 
for the KPZ equation are 2v — 1 and D — 1. Since the asymmetry equals /3, the 
WASEP height (4.12) should hence be compared with /3hp(x,t), see (2.11). Then, 
stretching somewhat the results in |9], for the WASEP limit it holds 

lim ^e(3h e ( [e^x] , e~ 2 t) + yhe' 1 - ±(3H - \og(2^e//3) = ph p {x, t) . (4.13) 



£->0 



The BG result (3.4) corresponds to the particular case (3 = 1 and non-singular initial 
data. 

Combining (4.13) and (3.9) one arrives at an exact formula for the one-point 
distribution of hp which reads 

F(/3hp(x,t) + ±(3H + (x 2 /2t) - 21og/3 < 7t s) = F t (s) . 



(4.14) 



To generalize to arbitrary A, u, D one starts from the definition (2.10) for the narrow 
wedge solution h(x,t) of the KPZ equation and uses the scale invariance of white 
noise and of Brownian motion to obtain 

r t/2u 

(\/2u)h(x,t/2u) = logE^iexp [(\VD/2v) J r){b{2vs), (t/2v) - s)ds] :<?(&(*)) 

= log E x (: exp [(A ^D/(2vf 2 ) J* V (b(s),t - s)ds] :5(b(t)) 
= ah a (x, t) 

with a fixed to 



(4.15) 



a 



{2v)-' i / 2 \D 1 ' 2 . 



(4.16) 



Inserting in (4.14) one obtains the one-point distribution of the KPZ equation 



with narrow wedge initial data, 

F((\/2u)h(x,t/2u) + ^( 7t ) 3 + (x 2 /2t) - 2 log a < lt s) = F t (s) 



(4.17) 



where 7t = 2 1 / 3 (a 4 t) 1 / 3 . The identity (4.17) is the central result of our contribution. 



Remark: After posting this article it was brought to our attention that G. Amir, I. 
Corwin, and J. Quastel independently obtained the identity (4.17) and the formula 
(3.13) for the distribution function F t . 
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5 Some properties of the one-point distribution 



(i) Scale invariance. In terms of the original parameters we define a as in (4.16) 
and introduce dimensionless space, time, and height through 

X = a 2 x, T = 2uaH, H = (X/2u)h . (5.1) 

Then the coefficients in the KPZ equation become A = 1, v — |, and D — 1, i.e. 
H(X, T) satisfies 

In these variables ( |4.17 ) reads 

¥(H(X,T) + ±T + (X 2 /2T) - 2 log a < ^ T s) = F T (s) (5.3) 

with 7t = 2~ 1//3 T 1 / 3 . The distribution function Ft(s) depends on T through 
substituting j t by jt- The solution (4.17) satisfies the scale invariance of the 
KPZ equation, the term 2 log a resulting from the scaling of the initial condition 
Z(x,0) = 6(x), see ( pH| . 

(it) Stationarity. Since for our initial conditions the height profile is curved, no 
stationarity in x is expected. However, in the case of PNG droplet growth it was 
noted that height fluctuations become stationary once the systematic curvature is 
subtracted. The same feature appears for the KPZ equation. Using dimensionless 
units for simplicity, we set 

h(x,t) =logE (:exp [J rj(b(s), s)ds\: 6(b(t) - ac)) , (5.4) 

where, with no modification of the statistics, the Brownian motion b(s) has been 
time-reversed to b(t — s). Let us introduce the Brownian bridge w(s) running from 
at s = to x at s = t. In terms of Brownian motion, 

w (s) = b(s) - ^(b(t) -x) . (5.5) 

Denoting by the expectation over a Brownian bridge with endpoints 0,x, we 
define t 

h(x,t) = logE™(:exp [ / v (w(s), s)ds]:) , (5.6) 

Jo 

which satisfies the ^-independent normalization 

(e^) = 1 . (5.7) 

Using scale invariance, one checks that x i— > h(x,t) is stationary as a stochastic 
process for fixed t. On the level of the one-point distribution this property is reflected 
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by F t (s) being independent of x. 

(in) Ground state energy of the 5-Bose gas. As first shown by McGuire 16 , in units 



given by 



of (2.18) with 2v — 1 the ground state energy of the 5-Bose gas with n particles is 

(5.8) 

-tE n 



E n = -±n(n 2 -l)a 4 



By (2.18), 



(Z(0,i) n ) = (0\e- tHn \0) 



(5.9) 



for large t. Of course, there are subleading terms, which however will be ignored at 
this stage. From the exact solution one has 



Hence 



F(ah a (0,t) + ±aH - 2 log a < s) = F t (s/ lt ) 

E(exp [n(ah a (0,t) + ±aH - 2 log a)]) 

= J ds J due ns e s e~^ u exp [-e s e~^ u ]g t ( 

= n\ I due nitu g t {u) 



(5.10) 



(5.11) 



with gt(u) defined in (3.14) 



7t diverges as t 1 ' 3 for t — > oo. Hence the integral in (5.11) is sensitive only to 



the right tail of gt(u). Furthermore, by the argument before (3.15), gt(u) tends to 
a right tail as exp[— |u 3//2 ], see e.g. [17]. 

POO 

/ dwexp[— fw 3//2 + n^ t u] = exp[^aHn 3 
Jo 



F^ w (m), which has a right tail as exp[— |u 3 / 2 ], see e.g. [17]. Therefore to leading 
order 



(5.12) 



for t — > oo. Combining with (5.11) 



lim — logE(exp \nah a (0, t)] ) 

t^oo t 



hn(n 2 — l)a 
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(5.13) 



in agreement with (5.8), (5.9). 



As discussed in [18], Chapters 9 and 10, the connection to the 5-Bose gas was 
used to support the scaling exponent 1/3 and to study the right tail of the height 
distribution. We added here the computation of the exact prefactor 4/3 for the 
stretched exponential. Note that the left hand tail of the Tracy- Widom density 
decays as exp[— |it| 3 /12]. 

6 Numerical evaluations 



According to (4.17) the one-point height statistics is governed by the probability 
density 

d 



Pt{s) 



ds 



-FAs] 



ltd 



■yt(s-u) 



exp [ 



_ p 7t (s-u) 



]g t (u)du, 



(6.1) 
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where gt(u) is defined in (3.14). This expression is sufficiently concrete to allow for 
some exploratory numerical studies. In Appendix A we establish that gt{s) can be 
approximated for large t as 

9t(s) = Ptw(s) + 9t(s) + 0(^ 6 ) (6.2) 



with the first order correction gt(s) given in (A. 17). Correspondingly one has the 
first order approximation to Pt(s) as 

pS 1} ( S ) = J 7t e*<-«> exp [ - e^"")] (p TW (u) + &(u))du . (6.3) 

ljt can be expressed through quantities related to the Tracy-Widom distribution, 



see (A. 30), and tabulated in 20 . The evaluation is displayed in Fig. 1, wher e 
the vertical axis corresponds to 74 = 1. In Fig. 2 we plot as defined in (6.3). 
Pt turns negative at 7$ = 5, which signals the border of validity of the first order 
approximation. As shown in Fig. 3, at 7$ = 2 the difference Pt — pxw has two 
nodes while the difference p t — Ptw as displayed in Fig. 5 below has only a single 
node, again indicating qualitative differences. 

A more demanding approach is a numerical evaluation of the determinants defin- 
ing g t . For this purpose the kernels of B t and P\i are computed using Mathematica. 
One chooses a properly adjusted equally spaced grid of 120 points for each variable 
resulting in a matrix {Mjj} ij= i i ... i i2o- The step size for s is 0.1 covering the interval 
[—6, 6]. The projection P s means to restrict M to i,j = (60 + 10s), . . . , 120. The 
respective determinant is then evaluated for various s, t again using Mathematica. 
To reach higher precision and earlier times it would be advisable to use the more 



refined methods of 21 



Pt is displayed for 7 t = 2, 5, 10 in Fig. 4. The dotted line is the limiting Tracy- 
Widom density. For better visibility, in Fig. 5 we plot the difference Pt(s)— Ptw( s ) a t 
the same times. In Fig. 6 we test directly the precision of the first order correction 
by plotting g t at j t = 10 and g t — Ptw( s )- While there is qualitative agreement 
we would have expected a much smaller error, in particular at the left minimum. 
A much cruder comparison is to plot merely the first moments of pt and pi\ As 
displayed in Fig. 7 the agreement is fairly precise beyond 74 = 5. 

As argued in the Introduction, features of our exact solution should be visible 
also in the corner growth model with finer and finer details being revealed as q is 
decreased from 1 to | + (3y/e. To support this claim we carried out a few Monte 
Carlo simulations. Time steps have a sufficiently fine discretization and random 
sequential updating is used. We measure the total number of signed jumps between 
sites and 1. For q — 1 at t — 10 3 MC steps the height distribution has an effective 
range of 30 equally spaced points. Except for the shifted mean, the Tracy-Widom 
density is an accurate approximation [3]. For earlier times the height takes so few 
values that correction to pxw does not seem to be a meaningful notion. In Fig. 8 
we plot the corresponding result for q = 0.6 at t — 1024 MC steps with an average 
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over 10 4 realizations. We compare with the exact solution at 7 t = 0.94. Clearly, 
this is a better approximation than the Tracy- Widom distribution indicated by the 
dashed line. 
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A Appendix: Expansion in 1/t 

Expansion of the determinants. We follow as closely as possible the notation 
in 1 19 1, simply referred to as TW. As only difference we set 



ij)(x) = Ai(x) (A(x) in TW) (A.l) 

and 

ip x (x) = ij){x + A) , ipo(x) = V(a;) • (A.2) 

We define the operator 

B t = K + C t , (A.3) 

with the kernels 

poo 

K(x,y)= / d\M%)My) , (A.4) 



o 



formerly denoted by Kai, and 

POO 

C t (x,y)= / dA(e^ A -l)- 1 (^ A (x)^A(l/)-^-A(^-A(z/)), (A.5) 



and the unnormalized projection with kernel ip(x)ip{y). 
We consider the difference 

g t (s) = det(l -K-Ct + P*)- det(l — K — C t ) . (A.6) 



The projection P s in (3.14) is taken into account by defining the determinant in 
L 2 ([s,oo)) with scalar product (•,•). For fixed s we want to expand g t (s) in 7^, 
which means to regard Ct as a small perturbation. 
One starts from 

g t = det(l - K - CM, (1 - K - a)'V) 

= det(l - K) det(l - (1 - K)- x C t )($, (1 - K - C t )~V> , (A.7) 
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expands the Fredholm determinant as 

det(l - (1 - K)- l C t ) = 1 - tr((l - KY l C t ) + 0{C 2 t ) , (A.8) 

and writes the resolvent as 

(i-K-Ct)- 1 ^) = (V, (l-K)-^)+{iP, {l-KY^til-KY^+OiCD . (A.9) 

The leading term equals F'/F, see [12], where 

F = F TW (s) = det(l - P s K Ai P s ) = det(l - AT) . (A. 10) 

The term linear in is 

& = F{(if>, (1 - A) _1 C t (l - ATV> - (1 - ATV>tr((l - A) _1 C t )} 
= ^ - F' + 0(C f 2 ) . (A. 11) 



Inserting the definition of Ct from (A.5), one arrives at 



POO 

g t = F dX(e^ x - 1) _1 (G(A) - G(-X)) (A.12) 
J o 



with 



G(A) = (V, (1 - A)-Va)(^a, (1 - A)-V) - (ipx, (1 - A)-Va)(^, (1 - A)-V) • 

(A.13) 

Rescaling by jt one arrives at 

POO 

at = F^ 1 / dA(e A - \y\G{\/ lt ) - G{-X/ lt )) . (A.14) 
Jo 

Thus the expansion in 1/t amounts to a Taylor expansion of G(X) = G(X) — G(—X) 
at A = 0. 

One has 5(0) = 0, G"(0) = 0, since G is odd. Furthermore G'(0) = and 

G"'(0) = 6(«, (1 - ATVX^', (1 - if)-V> 

-W, (1 - A)-V')(^, (1 - iO~V» • (A.15) 

Hence, using that 

/■oo 

/ A 3 (e A -l)- 1 dA = 7r 4 /15, (A.16) 

it holds 

£ = 2 74 - V/15)F(<^, (i _ Ky l ip) (V', (1 - A)" V> 

-(v", (i - Aryxv>, (i - #rv» ■ (A.i7) 
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Since g t is of order j t 4 , in principle, the terms of 0{Cf) could yield a more 
slowly decaying contribution. By expanding to second order, one checks that this 



contribution is at least of order j t 5 . Thus (A. 17) is the leading order in the 1/t 
expansion. 

Link to the Tracy- Widom distribution. For a numerical computation of g t it 
is convenient to use the well tabulated F |20|. One introduces, see TW, 

Q=(l-K)- x i/>, P=(l-K)- 1 tf, 
q(s) = Q(s) , p(s) = P(s) , u(s) = ty,Q), v(s) = ty,P). (A.18) 



To compute (A. 17) one still needs {i/j",Q) and (ifj",P). 

(i) Computation of (ip",P). We claim that 

(^",P> = > 2 -p 2 ). (A.19) 

Proof: One starts from TW, below (2.14), 

P') = (tf, (1 - K)- l 4,») + [D, (1 - K)-^') . (A.20) 

On the other hand 

(^P'> = -^( S )P( S )-(^',P>. (A.21) 

Hence 

- 2«, (1 - Jf)" V) = ij'(s)p + {</,', [D, (1 - JO - V> • (A.22) 
From TW (2.13) and the definitions of R and p there, 

(^', [D, (1 - K)- 1 )^') = -v 2 + W,R)(rl/,p) (A.23) 

and from TW above (2.18), 

<V/, R) = -(i>'(s) - p) , <?//, P )=p, (A.24) 

which establishes the claim. □ 

(ii) Computation of (i/;",Q). We claim that 

p(s') 2 ds' + vu - qp . (A.25) 
Proof: Using TW (2.14) one starts from 

W,qr) = -i//{s)q-W,Q) 

= (if/,P)-(i//,Q)u + q(R,i(J) 

= (il/,P)-vu+q(p-if/(s)). (A.26) 
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Hence 

W,Q)=-W,P)+vu-qp. (A.27) 

Next we consider 

j s ^\P) = -^\s)p-W 1 R)p 

= -^'{s)p-{-^{s)+p)p=-p 2 , (A.28) 

which establishes the claim. □ 

Using the identities from TW, it holds 

u = F'/F, v = F"/2F, q 2 =u 2 -2v, p = q + qu . (A.29) 



In combination with (A. 17), (A. 18), (A. 19), and (A. 25), one arrives at 



g t {s) = (7rVl5) Tt - 4 F( S ) 

x(-2v{s)(q{s)p(s)+ J p(s') 2 ds'^ +u(s)(v(s) 2 +p(s) 2 )y (A.30) 

which is the starting point for numerical evaluations. In our notation, |20| tabulates 



F,q } q' } u. From (A.29) we then deduce p and v and only the computation of the 
integral remains. 
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Figure 5: The difference p t (s) — Ptw(s) for j t = 2 (o), 5 (+), 10 (•) based on the 
evaluation of the determinants. 
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Figure 7: First moment of pt 



and p\ 
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Figure 8: The exact KPZ density at 7 t = 0.94 (— ) and the PASEP Monte Carlo at 
q = 0.6, t = 1024 MC steps (•). The dashed line is the Tracy- Widom density. 



23 



